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Abstract. We survey some known facts and open questions concerning 
the global properties of 3+1 dimensional spacetimes containing a com- 
pact Cauchy surface. We consider spacetimes with an ^-dimensional Lie 
algebra of space-like Killing fields. For each I < 3, we give some basic 
results and conjectures on global existence and cosmic censorship. 
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1. Introduction 

In this review, we will describe some results and conjectures about the 
global structure of solutions to the Einstein equations in 3+1 dimensions. 
We consider spacetimes (M, g) with an ^-dimensional Lie algebra of space- 
like Killing fields. We may say that such spacetimes have a (local) isometry 
group G of dimension I with the action of G generated by space-like Killing 
fields. 

For each value I < 3 of the dimension of the isometry group, we state the 
reduced field equations as well as attempt to give an overview of the most 
important results and conjectures. We will concentrate on the vacuum case. 

In section we present the Einstein equations and give the 3+1 decom- 
position into constraint and evolution equations, cf. subsection 12.21 Due to 
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the gauge freedom in the Einstein equations, questions on the global prop- 
erties of solutions to the Einstein equations must be posed carefully. We 
introduce the notions of vacuum extension and maximal Cauchy develop- 
ment and state the uniqueness theorem of Choquet-Bruhat and Geroch, for 
maximal vacuum Cauchy developments. We also collect here a few basic 
facts about Killing fields on globally hyperbolic spacetimes. In subsection 
12.11 a version of the cosmic censorship conjecture appropriate for vacuum 
spacetimes, with compact Cauchy surface, is stated and in subsection 12.31 
we discuss a few gauge conditions which may be of use for the global exis- 
tence problem for the Einstein equations. Section [2] is ended by a discussion 
of a few PDE aspects of the Einstein equations which are relevant for the 
topic of this paper, cf. subsection 12.41 

In the cases I = 3 (Bianchi, cf. section |SJ) and a special case of I = 
2 (polarized Gowdy, cf. section |IJ), the global behavior of the Einstein 
equations is well understood, both with regard to global existence of the 
evolution equations and the cosmic censorship problem. For the general 
I = 2 case (local U(l) xU(l) symmetry, cf. section |IJ), there are only partial 
results on the global existence problem and the cosmic censorship problem 
remains open, although conjectures supported by numerical evidence give a 
good indication of what the correct picture is. 

In the cases I = 1 (U(l) symmetry, cf. section [SJ and 1 = (no symmetry, 
i.e. full 3+1 dimensional Einstein equations, cf. section H3) , the large data 
global existence and cosmic censorship problems are open. In the U(l) 
case conjectures supported by numerical evidence give a good idea of the 
generic behavior, and there is a small data semi-global existence result for 
the expanding direction due to Choquet-Bruhat and Moncrief |42[ I39j . 

For 3+1 Einstein gravity without symmetries the only global existence 
results known are the theorem on nonlinear stability of Minkowski space 
of Christodoulou and Klainerman, the semi-global existence theorem of 
Friedrich for the hyperboloidal initial value proble and the semi-global ex- 
istence theorem for spatially compact spacetimes with Cauchy surface of 
hyperbolic type, due to Andersson and Moncrief. These are all small data 
results, see section El for discussion. 

Due to the high degree of complexity of the numerical solution of the 
Einstein equations in 3+1 dimensions it is too early to draw any conclu- 
sions relevant to the asymptotic behavior at the singularity for the full 3+1 
dimensional Einstein equations, from the numerical studies being performed. 
However, an attractive scenario is given by the so called BKL picture, cf. 
section [7] for some remarks and references. 

The Einstein equations are derived from a variational principle, and can 
be formulated as a Hamiltonian system (or time-dependent Hamiltonian sys- 
tem, depending on the gauge), and therefore the Hamiltonian aspect of the 
dynamics should not be ignored, see eg. the work by Fischer and Moncrief 
on the Hamiltonian reduction of the Einstein equations, |67| and references 
therein. In fact, the Hamiltonian point of view on the Einstein equations has 
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played a vital role as motivation and guide in the development of the results 
discussed here. The notion that the Einstein evolution equations in terms of 
canonical variables can be viewed as the geodesic spray for a metric on the 
phase space (deWitt metric) modified by a curvature potential, is natural 
from the Hamiltonian point of view, and this picture has been relevant to 
the development of ideas on asymptotic velocity dominance, see sections 
and El 

In this review, however, we will concentrate exclusively on the differential 
geometric and analytical point of view. Even with this restriction, many 
important topics are left out and we make no claim of complete coverage. 
See also ^39 and jllUj for related surveys. 

Acknowledgements: I am grateful to Vince Moncrief for numerous con- 
versations on the topics covered here and for detailed comments on an early 
version. Thanks are due to Hakan Andreasson, Piotr Chrusciel, Jim Isen- 
berg and Alan Rendall and others for helpful comments. I am happy to 
acknowledge the hospitality and support of the Institute of Theoretical 
Physics, UCSB, and Institut des Hautes Etudes Scientifiques, Bures sur 
Yvette, where where part of the writing was done. 

2. The Einstein equations 

Let (M, g) be a smooth 4-dimensional Lorentz manifold 1 of signature 
— h++. The Lorentzian metric g defines a causal structure on M. For 
the convenience of the reader we give a quick review of the basic causality 
concepts in appendix lAl see (1651 Chapter 8]. |241 I5il 1132] for details. 

We will here consider only the vacuum case, i.e. the case when g is Ricci 
flat, 

Rab = 0. (2.1) 

Let M C M be a space-like hypersurface, i.e. a hypersurface with time-like 
normal T. We let ej be a frame on M and use indices k for the frame 
components. Let g, k be the induced metric and second fundamental form of 
M C M, where kij = (Vjej, T). Let t be a time function on a neighborhood 
of M. Then we can introduce local coordinates (t,x % ,i = 1,2,3) on M so 
that x 1 are coordinates on the level sets Mt of t. This defines the coordinate 
vector field dt of t. Alternatively we can let M t = i(i, M) where i:lxM-> 
M is a 1-parameter family of imbeddings of an abstract 3-manifold M. Then 
dt = i*d/dt where d/dt is the coordinate derivative on R. 

Define the lapse iV and shift X w.r.t. t by dt = NT + X. Assume that 
T is future oriented so that iV > 0. A 3+1 split of equation (|2.1j) gives the 



*We denote the covariant derivative and curvature tensors associated to (M, g) by 
V, Rabcd etc. All manifolds are assumed to be Hausdorff, second countable and C°° , and 
all fields are assumed to be C°° unless otherwise stated. 
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Einstein vacuum constraint equations 

R- \k\ 2 + {tikf = 0, (2.2a) 

Vitrk - V j kij = 0, (2.2b) 
and the Einstein vacuum evolution equations 

£>d t 9ij = -2Nkij + Cxgij, (2.3a) 
C dt k i:j = -ViVjN + N{R i:j + trkhj - 2ki m k n ]) + Cxhj. (2.3b) 

where £,g t denotes Lie-derivative 2 w.r.t. dt- In case [c?t,ei] = 0, Cg t can be 
replaced by dt- 

A triple (M, g, k) consisting of a 3-manifold M, a Riemannian metric g 
on M and a symmetric covariant 2-tensor A; is a vacuum data set for the 
Einstein equations if it solves (|2.2|) , 

Definition 2.1. Let (M,g,k) be a vacuum data set. 

(1) A vacuum spacetime (M, g) is called a vacuum extension of (M, g, k) 
if there there is an imbedding i with time-like normal T of (M, g, k) 
into (M,g) so that g = i*g and k = — i*(VT). 

(2) A globally hyperbolic vacuum spacetime (M,g) is called a vacuum 
Cauchy development of (M, g, k) if there is an imbedding i with 
time-like normal T of (M, g, k) into (M, g) so that i(M) is a Cauchy 
surface in (M,g), g = i*g and k = — i*(VT). If (M,g) is maximal 
in the class of vacuum Cauchy developments of (M,g, k) then (M,g) 
is called the maximal vacuum Cauchy development (MVCD) 
of (M,g,k). In the following, when convenient, we will identify M 
with i(M). 

The Einstein vacuum equations are not hyperbolic in any standard sense 
due to the coordinate invariance ("general covariance") of the equation 
R a b = 0. Nevertheless, the Cauchy problem for the Einstein vacuum equa- 
tion is well posed in the following sense. 

Theorem 2.2 (Choquet-Bruhat and Geroch i40 ]). Let (M,g,k) be a vac- 
uum data set. Then there is a unique, up to isometry, maximal vacuum 
Cauchy development (MVCD) of (M,g,k). If (ft : M — > M is a diffeomor- 
phism, the MVCD of(M,(p*g,(j)*k) is isometric to the MVCD of(M,g,k). 

The proof relies on the fact that in spacetime harmonic coordinates, 
OgX a = 0, the Ricci tensor is of the form 

5%j> = -\u- g g ap + S aP [g, dg], (2.4) 

where \Dg is the scalar wave operator in (M,g). Hence the Einstein vacuum 
equations in spacetime harmonic coordinates is a quasi-linear hyperbolic 

2 for a tensor b on M, with bij = &(ej,e,), we have Ld t bij = dtipij) — b([dt, et], ej) — 
b(ei, [d u ej]) 
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system and therefore the Cauchy problem 3 for R^ = is well posed and 
standard results give local existence. One proves that if the constraints and 
gauge conditions are satisfied initially, they are preserved by the evolution. 
This together with a Zorn's lemma argument gives the existence of a MVCD. 
Uniqueness is proved using the field equations to get a contradiction to the 
Haussdorff property, given a pair of non-isometric vacuum Cauchy develop- 
ments, which are both maximal w.r.t. the natural partial ordering on the 
class of Cauchy developments. 

A spacetime (M,g) is said to satisfy the time— like convergence con- 
dition (or strong energy condition) if 

RabV a V b > 0, for all V with g ah V a V h < 0. (2.5) 

Globally hyperbolic spacetimes with compact Cauchy surface and satisfying 
the time-like convergence condition are often called "cosmological space- 
times" in the literature, following [23]. Here we will use the term spatially 
compact to refer to the existence of a compact Cauchy surface. A spacelike 
hypersurface (M,g) in (M,g) has constant mean curvature if Vitr g k = 0, 
cf. subsection 12 . 31 below . 

We end this subsection with by stating a few facts about Killing fields. 

Proposition 2.3 (|68j). Let M be a compact manifold and let (M,g,k) be 
a constant mean curvature vacuum data set on M with MVCD (M,g). Let 
Y be a Killing field on (M,g) and let Y = Y±T + Ki be the splitting ofY 
into its perpendicular and tangential parts at M. Then (Y±,Yn) satisfy the 
conditions 

(1) Y± = 0, £>Yn9 = 0, Cy,,k = 0, in case g is non-flat or k ^ 0. 

(2) Y± is constant and £>y«9 = if g is flat and k = 0. 

On the other hand, given Yx,Y\\ on M satisfying conditions'® above, there 
is a unique Killing field Y on M , with Y = Y±T + Yu on M . □ 

Proposition 2.4. Let (M,g) be a globally hyperbolic spacetime. 

(1) Assume that (M,g) satisfies the time-like convergence condition and 
contains a compact Cauchy surface M with constant mean curvature. 
Then either (M, g) is a metric product M xK or any Killing field Y 
on (M,g) is tangent to M. In particular, if (M,g) is vacuum and 
has a nonzero Killing not tangent to M , then (M, g) is flat. 

(2) Assume a compact group G acts by isometries on (M,g). Then the 
action of G is generated by space-like Killing fields and (M, g) is 
foliated by Cauchy surfaces invariant under the action of G. 

(3) Assume that (M, g) is 3+1 dimensional. Let M be a Cauchy surface 
in M , let Y be a Killing field on M and assume Y is strictly spacelike, 
g(Y,Y) > 0, on M. Then Y is strictly spacelike on M. 



Note that the Einstein equations in spacetime harmonic gauge should be viewed as 
an evolution equation for (g, k, N, X). 
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Proof. Point ^ is a well known consequence of the uniqueness result for 
constant mean curvature hypersurfaces of Brill and Flaherty [35,, cf. |124j . 
Point 121 is essentially |27| Lemma 1.1]. For the proof, note that as G is 
compact we can construct a G invariant time function on M by averaging 
any global time function t on M w.r.t. the G action, cf. the proof of |27[ 
Lemma 1.1]. The level sets of the averaged time function are Cauchy surfaces 
and are invariant under the action of G. The result follows. 

The following argument for point is due to Alan Rendall 4 . Let N = 
{p E M : g(Y,Y) = 0} and assume for a contradiction N is nonempty. 
Choose p £ N and a time function t on M so that t(M) = and t(p) > 0. 
Let A denote the intersection of the past of p with the future of M and let 
t\ = inf {£(</) : q £ AnN}. The set INT is closed and by global hyperbolicity A 
is compact and hence t\ > and there is a q £ An'N with t(q) = t\. If Y(q) 
is nonzero and null, then using the equation Y a V a (Y fe Y,) = which holds 
since Y is Killing, gives a null curve of points in N where Y is null. Following 
this into the past, shows that there is a q G N with < fi, which gives 
a contradiction. In case Y{q) = 0, the linearization of Y acts by isometries 
on T q M, and as the sphere of null directions in T q M is two dimensional 
it leaves a null direction fixed. Using the exponential map shows that the 
action of Y near q leaves a null geodesic invariant, along which Y must be 
zero or null. This leads to a contradiction as above. □ 

2.1. Cosmic Censorship. Theorem 12.21 proves uniqueness of the MVCD 
of a given data set (M,g, k). However, examples show that the MVCD may 
fail to be maximal in the class of all vacuum extensions, i.e. there exist 
examples of vacuum data sets (M, g, k) with vacuum extensions (M, g) such 
that the MVCD of (M, g, k) is a strict subset of (M,g). 

Example 2.5. Consider the n + 1-dimensional Minkowski space l"' 1 with 
metric rj = —dt 2 + (dx 1 ) 2 + • • • + (dx n ) 2 , let I + ({0}) be the interior of the 
future light cone. 7 + ({0}) is globally hyperbolic with the hyperboloids as 
Cauchy surfaces, and with the mantle of the light cone as Cauchy horizon. 
Let r be a cocompact discrete subgroup of the Lorentz group SO(n, 1). Then 
the quotient space M = r\/ + ({0}) is a globally hyperbolic, spatially compact 
spacetime. By choosing p to be the Lorentzian distance from the origin, 
we get g = —dp 2 + p 2 ^ where 7 is the standard hyperbolic metric on the 
compact quotient M = T\H n . (M,g) is the MVCD of the vacuum data set 
(M, 7 ,- 7 ). 

In case n = 1, H 1 = R, and a fundamental domain for V can be found 
which intersects the null boundary o// + ({0}) in an open interval. Therefore 
if n = 1, there is a nontrivial extension of M, which is still flat, but which 
fails to be globally hyperbolic, cf. figure^ This spacetime is known as the 
Misner universe. The maximal extension is unique in this case. 



private communication, 1999 
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Figure 1. The 1+1 dimensional Misner universe, showing 
a few orbits of S0(1, 1) (dashed lines) and orbits of a dis- 
crete subgroup T C S0(1, 1) (circles). The boundary of a 
fundamental domain for T is plotted with solid lines. 

If n > 1, the ergodicity of the geodesic flow on (M, 7) can be shown to 
prevent the existence of a spacetime extending (M, g) |lf)4j . 

Higher dimensional examples of flat globally hyperbolic spatially compact 
spacetimes which admit nontrivial non-globally hyperbolic extensions can be 
constructed by taking products of the n = 1 Misner universe with the flat 
torus. 

A maximal vacuum extension may be non-unique, as is shown by the 
Taub-NUT example, cf. 53 j. If the MVCD of a vacuum data set is not 
a maximal vacuum extension, any extension of it must fail to satisfy the 
intuitively reasonable causality requirement of global hyperbolicity. 

According to physical intuition, causality violations should be rare. This 
leads to the idea of cosmic censorship, essentially due to Penrose, see |133| 
for discussion. One way of stating this, relevant to the class of spacetimes we 
are concerned with here, is the following form of the strong cosmic censorship 
conjecture. 

Conjecture 1 (Strong Cosmic Censorship). Let M be a compact manifold 
of dimension 3. Then for generic vacuum data sets (M,g,k), the maximal 
vacuum Cauchy development of (M, g, k) is equal to the maximal vacuum 
extension of (M,g, k). 

In the case of asymptotically flat spacetimes (describing isolated systems 
in general relativity), the so-called weak cosmic censorship conjecture states 
that naked singularities (i.e. singularities which can be seen by an observer 
at infinity) should not occur generically, see the review paper by Wald |166j 
for a discussion of the status of the weak cosmic censorship conjecture. The 
work of Christodoulou, see and references therein, see also the discussion 
in |166[ §5], establishes weak cosmic censorship in the class of spherically 
symmetric Einstein-scalar field spacetimes, but also gives examples of initial 



8 L. ANDERSSON 

data such that the Cauchy development has a naked singularity. For earlier 
surveys on the strong cosmic censorship conjecture, see [25j and |l9] . 

The Penrose inequality, giving a lower bound on the ADM mass in terms 
of the area of a horizon in black hole spacetimes, was derived by a heuristic 
argument assuming the validity of the weak cosmic censorship conjecture. 
The proof of the Riemannian version of the Penrose inequality by Huisken 
and Ilmanen |91| 192). gives indirect support for the conjecture. 

Let (M, g) be a spacetime and let M C M be a space-like hypersurface. 
The Cauchy horizon H{M) is the boundary of the domain of dependence 
D(M), cf. Appendix for definition. If M is compact without boundary, 
then every point of H + {M) lies on a past inextendible null geodesic and 
every point of H~{M) lies on a future inextendible null geodesic, where 
H + (M) and H~{M) are the future and past components of H(M), respec- 
tively. 

Let (M, g) be a maximal vacuum extension of a vacuum data set (M, g, k) 
with M compact, and let D(M) C M be the MVCD of (M, g, k). If D(M) ^ 
M, then the Cauchy horizon H(M) is nonempty. One approach to SCC is to 
study the geometry of Cauchy horizons in vacuum spacetimes and to prove 
rigidity theorems as a consequence of extendibility of D(M). 

Isenberg and Moncrief proved for analytic vacuum or electrovac space- 
times, with analytic Cauchy horizon H(M), that under the additional as- 
sumption that H(M) is ruled by closed null geodesies, there is a nontrivial 
Killing field which extends to D(M), see jlUUI 1130] . This result was gener- 
alized to the C°° case by Friedrich et. al. As spacetimes with Killing 
fields are non-generic, this may be viewed as supporting evidence for the 
SCC. 

In the class of Bianchi spacetimes (i.e. spatially locally homogenous space- 
times, cf. section EJ), it has been proved by Chrusciel and Rendall [55] . 
generalizing work by Siklos |147j in the analytic case, that any C°° Bianchi 
spacetime which contains a compact locally homogenous Cauchy horizon is 
a Taub spacetime, cf. section 03 1)3.7)1 for definition. This result may be 
viewed as a version of SCC in the class of Bianchi spacetimes. In this con- 
text, it is worth mentioning that work by Chrusciel and Galloway 52 gives 
examples which indicate that Cauchy horizons may be non-differentiable, 
generically. 

Conjecture 2 (Bartnik [231 Conjecture 2]). Let (M,g) be a spatially com- 
pact globally hyperbolic spacetime satisfying the time-like convergence con- 
dition 12. fy) . Then if(M,g) is time-like geodesically complete, (M,g) splits 
isometrically as a product (R x M, —dt 2 + g). 

If the Bartnik conjecture [2] is true, then any vacuum, globally hyperbolic, 
spatially compact spacetime, is either flat and covered by R x T 3 or is has 
an inextendible time-like geodesic which ends after a finite proper time, i.e. 
it is time-like geodesically incomplete. A sequence of points approaching 
the "end" of a finite length inextendible geodesic is often thought of as 
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approaching a singularity. See [20] for a discussion of the status of the 
Bartnik conjecture. 

Inextendibility of D(M) can be detected by monitoring the asymptotic 
behavior of curvature invariants such as the Kretschmann scalar k, defined 
by k = R a b C( iR abcd ■ If k blows up along causal geodesies, then D(M) fails to 
be extendible, and therefore proving blowup for k for generic spacetimes is 
an approach to proving SCC. This is the method used in the proof of cosmic 
censorship for the class of polarized Gowdy spacetimes [5lJ, cf. section 01 
and is likely to be important also in the cases with less symmetry. The 
structure of the horizon and extensions in the polarized Gowdy class can 
be very complicated as shown by the work of Chrusciel et. al., see |49| for 
discussion, see also [SS]. It was proved by Ringstrom, cf. Theorem 13. 11 that 
for vacuum Bianchi spacetimes of class A, either the spacetime is Taub, cf. 
section |31 or k blows up at the singularity. 

2.2. The evolution equations. A solution to the vacuum Einstein evolu- 
tion equations with initial data is a curve t (g(t), k(t),N(t),X(t)) defined 
on some interval (To,Ti), satisfying (j2.3j) . 

Every regular solution (g, k, N, X) to the vacuum evolution equations 
(|2.3j) with initial data solving the vacuum constraint equations, gives a vac- 
uum spacetime. This is due to the fact that the constraint quantities 

B = R+(trk) 2 - \k\ 2 
Di = Vitik - V j kij 

evolve according to a symmetric hyperbolic system and energy estimates 
together with an application of the Gronwall inequality allow one to show 
the the constraints are satisfied during the time of existence of the solution 
curve. Now the fact that the Einstein vacuum equation is equivalent to 
the system of constraint and evolution equations shows that the spacetime 
(M, g) constructed from the curve (g, k, N, X) by letting M = (To, T\) x M, 
and setting 

g = -N 2 dt 2 + g ij (dx i + X i dt)(dx j + X j dt) 

is a solution to the Einstein vacuum equations (|2.1|) . 

Note that in order for the solution to be well defined, it is necessary to 
specify the lapse and shift (N,X), either as functions on spacetime M x 
(Tq,Tx) or as functions of the data, N = N[g,k],X = X[g, k]; this may be 
viewed as a gauge fixing for the Einstein equations. 

The choice of lapse and shift is crucial for the behavior of the solution 
curve. In particular, a foliation constructed for a particular choice of N, X 
may develop singularities which are not caused by any singular or irreg- 
ular nature of the Cauchy development. Consider for example the Gauss 
foliation condition N = 1,X = 0. Then the hypersurface M flows in the 
direction of its unit normal and Mt is simply the level set of the Lorentzian 
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distance function t(p) = d(M,p). The foliation {Mt} will develop singulari- 
ties precisely at the focal set of M, which in general will be nonempty, even 
in Minkowski space. 

Many authors have considered hyperbolic reformulations of the Einstein 
equations, see the paper by Friedrich j2J for discussion, see also [231131113] 
for related work. The development of singularities for hyperbolic systems 
presents a serious obstacle to the numerical treatment of the Einstein evolu- 
tion equations using hyperbolic reformulations, see p] |2j for discussion and 
examples. It is therefore necessary to consider also gauges which make (|2.<3|) 
into an elliptic-hyperbolic system. 

2.3. Constant mean curvature foliations. A particularly interesting 
choice of gauge condition for the lapse function is given by the constant 
mean curvature (CMC) condition 

VjtrA; = 0, 

i.e. the level sets Mt of the time function t are assumed to be hypersurfaces of 
constant mean curvature in (M, g). If (M, g) is globally hyperbolic, spatially 
compact and satisfies the time-like convergence condition (|2.5|) . then for 
r S R, either there is at most one Cauchy surface with mean curvature 
trfc = r or (M, g) splits product, cf. [SSI- 

This indicates that the mean curvature trk may be useful as a time func- 
tion on (M, g) in the spatially compact case. Setting t = trk leads, using 
(|2.2H2.3|) to the lapse equation 

-AN +\k\ 2 N = 1. (2.6) 

Local wellposedness for the Einstein evolution equations in CMC and more 
general elliptic time gauges given by a lapse equation of the form — AN + 
|/c| 2 ./V = NT(h) for given spacetime functions h has been proved by Choquet- 
Bruhat and York 03] . We may call this gauge the prescribed mean curvature 
(PMC) gauge. The proof is based on writing a wave equation for the second 
fundamental form, a technique which was also used in the work of Choquet- 
Bruhat and Ruggeri [13! on harmonic time gauge. The prescribed mean 
curvature condition, has been used in the U(l) x U(l) symmetric case by 
Henkel [HE] who proved global existence in this gauge. An important advan- 
tage of using PMC gauge is that one avoids the problem of proving existence 
of a CMC Cauchy surface. 

Equation 1)2.6(1 supplemented by an elliptic shift gauge makes the Ein- 
stein evolution equations ()2.3() into an elliptic-hyperbolic system of evolution 
equations. Let g be a given smooth metric on M, and let V 1 = g mn (T' l rnn — 
T l mn )- Then — V 1 is the tension field of the identity map Id : (M, g) — > (M, g), 
and the CMC condition coupled with the spatially harmonic coordinate con- 
dition 



V i = 
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gives a new elliptic gauge condition for the Einstein equations, CMCSH 
gauge ^0]- In jlTTj . it is shown that the Einstein equations in CMCSH gauge 
forms an elliptic-hyperbolic system for the Einstein equations, which is well 
posed in H s x H s ~ x , s > n/2 + 1. 

The maximal slicing condition trk = is of interest mainly for the asymp- 
totically flat case. This was used in the proof of the nonlinear stability of 
Minkowski space by Christodoulou and Klainerman [Ha]) cf. the discussion 
in section El Due to the "collapse of the lapse" phenomenon, see , the 
maximal foliation is not expected to cover the whole MVCD except in the 
small data case. See |124| for a discussion of maximal slices. Asymptoti- 
cally flat spacetimes satisfying certain restrictions on the causal structure 
are known to contain maximal hypersurfaces The existence of CMC 

hypersurfaces in asymptotically flat spacetimes was considered in [SJ. 

The mean curvature operator satisfies a geometric maximum principle, see 
7 for a proof of this under weak regularity. This allows one to use barriers 
to prove existence of constant mean curvature hypersurfaces. A spacetime 
is said to have crushing singularities if there are sequences of Cauchy 
surfaces with mean curvature trk tending uniformly to ±oo. Gerhardt |78j 
proved, using a barrier argument, that any spacetime satisfying (|2.5|) with 
crushing singularities is globally foliated by CMC hypersurfaces. These facts 
indicate that the CMC foliation condition is an interesting time gauge for 
the Einstein evolution equation. See also [7^1 ESI HH H31 EH! f° r results 
relevant to existence of CMC hypersurfaces. 

Let R[g] be the scalar curvature. A 3-manifold M is said to be of Yamabe 
type —1 if it admits no metric with R[g] = (and hence no metric with 
nonnegative scalar curvature), of Yamabe type if it admits a metric with 
R[g] = but no metric with R[g] = 1 and of Yamabe type +1 if it admits a 
metric with R{g] = 1, cf. 67, Definition 9]. 

If the Cauchy surface M is of Yamabe type —1, it follows from the con- 
straint equation that (M, g) cannot contain a maximal (i.e. trk = 0) Cauchy 
surface and therefore one expects that (if the dominant energy condition 
holds) the maximal time interval of existence for (|2,3j) in CMC time is of 
the form (after a time orientation) (— oo,0) with r /* corresponding to 
infinite expansion. 

If M is of Yamabe type 0, then one expects that either the maximal 
CMC time interval is (— oo,0) (possibly after a change of time orientation) 
or (M, g) splits as a product, and therefore in the vacuum case is covered by 
IxT 3 with the flat metric. Finally in case M is of Yamabe type +1, one 
expects the maximal CMC time interval to be (— oo,oo), i.e. the spacetime 
evolves from a "big bang" to a "big crunch". This is formalized in the 
"closed universe recollapse conjecture" of Barrow, Galloway and Tipler [2U]. 

Conjecture 3 (Constant mean curvature foliations). Let M be a compact 
3-manifold and let (M,g, k) be a vacuum data set on M , with constant mean 
curvature. The Cauchy problem for the Einstein vacuum evolution equations 
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with data (M, g, k) has global existence in the constant mean curvature time 
gauge, i.e. there is a CMC foliation in the MVCD (M,g) of (M,g,k), 
containing M , with mean curvature taking all values in (—00, 00) in case M 
has Yamabe type +1 and in case M has Yamabe type or —1, taking all 
values in (— 00, 0) (possibly after a change of time orientation). 

Remark 2.1. Conjecture^ has been stated in essentially this form by Ken- 
dall [135, Conjecture 1], see also Eardley and Moncrief ,1291 Conjecture C2] 
for a closely related statement. 

Note that as r \ —00, the past focal distance of the (unique) CMC surface 
with mean curvature r tends to zero, and hence the foliation exhausts the 
past of M . It follows that in case M has Yamabe type +1, then if Conjecture 
[3] is true, (M, g) is globally foliated by CMC hyper surf aces. In case M has 
Yamabe type or —1 on the other hand, there is the possibility that the CMC 
foliation does not cover all of M , due to the fact that as the mean curvature 
t /* 0, the CMC hyper surf aces are expected to avoid black holes, by analogy 
with the behavior of CMC and maximal hypersurfaces in the Schwarzschild 
spacetime. See |135j for further remarks and conjectures related to this. 

If one were able to prove Conjecture El then as remarked in |129j . this 
would give the possibility of attacking the Cosmic Censorship Conjecture 
using PDE methods. There are no known counterexamples to Conjecture El 
for vacuum spacetimes. However, Isenberg and Rendall 103 give an exam- 
ple of dust spacetimes, not covered by a CMC foliation. Bartnik (23] gave 
an example of a spatially compact, globally hyperbolic spacetime satisfy- 
ing the time-like convergence condition, which contains no CMC Cauchy 
surface. It is an interesting open question whether or not similar counter 
examples are possible in the vacuum case. The gluing results for solutions 
of the constraint equations, see jHZl EE1 EE] m &y lead to the construction of 
such examples. 

The CMC conjecture El has been proved in a number of cases for space- 
times with symmetry, in the sense of the existence of a group G of isometries 
acting (locally) on (M, g) by isometries and with space-like orbits. In the 
case of Bianchi IX, cf. section El which has Yamabe type +1, the closed 
universe recollapse conjecture and consequently Conjecture El was proved 
by Lin and Wald |121j . In the case of 2+1-dimensional vacuum spacetimes 
with cosmological constant, the conclusion of Conjecture El is valid |11| . In 
the 2+1 case, the CMC foliations are global. 

The work of Rendall and Burnett, see |36| and references therein, proves 
under certain restrictions on the matter that a maximal, globally hyperbolic, 
spherically symmetric spacetime, which contains a CMC Cauchy surface 
diffeomorphic to S 2 x S , is globally foliated by CMC hypersurfaces with 
mean curvature taking on all real values. 

We end this subsection by mentioning the harmonic time gauge condition, 
D s t = or 

g ab r ° ab = 
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In case X = 0, this is equivalent to the condition TV" = \/ det g / \/ det e 
where e%j is some fixed Riemannian metric on M. The Einstein evolution 
equations with X = were proved to be hyperbolic with this time gauge 
by Choquet-Bruhat and Ruggeri |33]- This time gauge also appears in the 
work of Gowdy and is used in the analysis of the Gowdy spacetimes as well 
as in the numerical work of Berger, Moncrief et. al. on Gowdy and U(l) 
spacetimes, cf. sections H] and |3 

The papers of Smarr and York [1481 1149j contain an interesting discus- 
sion of gauge conditions for the Einstein equations. See also section |5] for 
a discussion of the spatial harmonic coordinate gauge and the survey of 
Klainerman and Nicolo jllUj for further comments on gauges. 

2.4. The Einstein equations as a system of quasi linear PDE's. As 

mentioned above, the Einstein vacuum equations in spacetime harmonic 
coordinates form a quasi-linear hyperbolic system of the form 

~n- g g a/3 + S af3 [g,dg} = 0. (2.7) 

The system (J2ZJl 

is a quasi-linear wave equation, quadratic in the first 
order derivatives dg and with top order symbol depending only on the field 
g itself. Standard results show that the Cauchy problem is well posed in 
Sobolev spaces H s x H s , s > n/2 + 1. This was first proved for the 
Einstein equations by Hughes, Kato and Marsden 90 . It is also possible to 
prove this for elliptic-hyperbolic systems formed by the Einstein evolution 
equations together with the CMC-spatial harmonic coordinates gauge, see 
|1U) . see also section for a discussion of the spatial harmonic coordinates 
gauge. 

Recent work using harmonic analysis methods by Bahouri and Chemin 
|18j . Klainerman and Rodnianski |115U114Tlll3| and Smith and Tataru [151 j . 
see also |157j has pushed the regularity needed for systems of the above type 
on W 1 ' 1 to s > (n+ l)/2 for n > 3. In order to get well posedness for s lower 
than the values given above, it is likely one needs to exploit some form of 
the null condition, [1541 11521 1157j . The null condition for equations of the 
form D^u = F[u, du] on Minkowski space states roughly that the symbol 
of the non-linearity F cancels null vectors. For a discussion of the null 
condition on a curved space background, see |153| 115(1] . Counter examples 
to well posedness for quasi-linear wave equations with low regularity data 
are given by Lindblad, see [122 . 

The standard example of an equation which satisfies the null condition is 
□u = d a ud/3UT] a ^ where □ is the wave operator w.r.t. the Minkowski metric 
7] on W 1 ' 1 . This equation is well posed for data in H s x H s ~ l with s > n/2 
|116j . and has global existence for small data for n > 3. On the other 
hand, the equation Du = (dtu) 2 which does not satisfy the null condition 
can be shown to have a finite time of existence for small data in the 3+1 
dimensional case. 
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For quasi-linear wave equations which satisfy an appropriate form of the 
null condition 89 , global existence for small data is known in 3+1 dimen- 
sions. The Einstein equations, however, are not known to satisfy the null 
condition in any gauge. In particular, it can be seen that in spacetime har- 
monic coordinates, the Einstein equations do not satisfy the null condition. 
However, the analysis by Blanchet and Damour [Slj of the expansion of so- 
lutions of Einstein equations in perturbation series around Minkowski space 
indicates that the logarithmic terms in the gravitational field in spacetime 
harmonic coordinates, arising from the violation of the null condition, can 
be removed after a (nonlocal) gauge transformation to radiative coordinates 
where the coordinate change depends on the history of the field. A similar 
analysis can be done for the Yang-Mills (YM) equation in Lorentz gauge. It 
may further be argued that the small data, global existence proof of Christo- 
doulou and Klainerman for the Einstein equations exploits properties of the 
Einstein equations related to the null condition. In a recent paper, Lindblad 
and Rodnianski |123j have introduced a weak form of the null condition and 
shown that the Einstein equations satisfy this. 

Global existence is known for several of the classical field equations such 
as certain nonlinear Klein-Gordon (NLKG) equations and the YM equation 
on M 3,1 (proved by Eardley and Moncrief, |63l I64p . The proofs for NLKG 
and the proof of Eardley and Moncrief for YM use light cone estimates 
to get apriori L°° bounds. The proof of Eardley and Moncrief used the 
special properties of YM in the radial gauge. This method was also used in 
the global existence proof for YM on 3+1-dimensional, globally hyperbolic 
spacetimes by Chrusciel and Shatah [20]. Klainerman and Machedon jl()9j 
were able to prove that the YM equations on M 3,1 in Coloumb gauge satisfy 
a form of the null condition and are well posed in energy space H 1 x L 2 . 
They were then able to use the fact that the energy is conserved to prove 
global existence for YM. See also for an overview of these ideas and 
some related conjectures. 

An important open problem for the classical field equations is the global 
existence problem for the wave map equation (nonlinear u-model, hyperbolic 
harmonic map equation). This is an equation for a map M™' 1 — ► N, where 
N is some complete Riemannian manifold, 

a v u A + Yi c {u)d a u B d p u c r] aP = 0. (2.8) 

Here T is the Christoffel symbol on N. 

The wave map equation satisfies the null condition and hence we have 
small data global existence for n > 3. Small data global existence is also 
known for n = 2. Further, scaling arguments provide counterexamples to 
global existence for n > 3, whereas n = 2 is critical with respect to scaling. 
For n > 2, global existence for "large data" is known only for symmetric 
solutions, and in particular, the global existence problem for the wave map 
equation (|2.8|) is open for the case n = 2. For the case n = 1, global existence 
can be proved using energy estimates or light cone estimates. Tao proved 
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|156j that the wave map equation in 2+1 dimensions with the sphere as 
target is globally well-posed for small energy. The proof is based on a new 
local well-posedness in H 1 x L 2 and conservation of energy. This result has 
been generalized to the case of hyperbolic target by Krieger |119j . Tataru 
|158| has proved that the wave map equation in n+1 dimensions is locally 
well-posed in H n / 2 x H n / 2 ^ 1 for general target imbedded in M. m . (Here 
H s is the homogenous Sobolev space. See |158j for the precise definition of 
well-posedness used). See also (1461 1159 for surveys. 

The above discussion shows that the situation for the wave map equation 
is reminiscent of that for the Einstein equations, cf. sections HI |SJ In partic- 
ular, it is interesting to note that equations of the wave map type show up in 
the reduced vacuum Einstein equations for the Gowdy and U(l) problems. 

3. BlANCHI 

Let (M, g) be a 3+1 dimensional spacetime with 3-dimensional local isom- 
etry group G. Assume the action of G is generated by space-like Killing 
fields and that the orbits of G in the universal cover of M are 3-dimensional. 
This means there is a global foliation of M by space-like Cauchy surfaces 
M with locally homogeneous induced geometry. Such spacetimes are known 
as Bianchi spacetimes. The assumption of local homogeneity of the 3— 
dimensional Cauchy surfaces means that a classification of the universal 
cover is given by the classification of 3-dimensional Lie algebras. 

Let e a , a = 0, . . . , 3 be an ON frame on M, with eo = u, a unit time-like 
normal to the locally homogeneous Cauchy surfaces, let j^b ^ e the commu- 
tators of the frame, [e a , e&] = 7„ fc e c . Let the indices i,j, k, I run over 1, 2, 3. 
We may without loss of generality assume that [e , £j] = where {£?,}f =1 is 
a basis for the Lie algebra g of G. 

Choose a time function t so that i a it a = 1, i.e. the level sets of t coincide 
with the group orbits. Restricting to a level set M of t, the spatial part of the 
commutators 7^. are the structure constants of q. These can be decomposed 
into a constant symmetric matrix n kl and a vector a^, 

lij = eijin kl + ai<5$ - a$. 

We will briefly describe the classification used in the physics literature, cf. 
HS1, [Ml §1-5.1]. 

The Jacobi identity implies n y aj = and by choosing the frame {e{\ to 
diagonalize n 13 and so that e\ is proportional to the vector a^, we get 

n lj = diag(ni, n 2 ,n 3 ), a« = (a, 0,0). (3.1) 

The 3-dimensional Lie algebras are divided into two classes by the condition 
a = (class A) and a / (class B). The classes A and B correspond 
in mathematical terminology to the unimodular and non-unimodular Lie 
algebras. If 712723 / 0, let the scalar h be defined by 

a 2 = hn2n^. (3.2) 



16 



L. ANDERSSON 



Class A 



Table 1. Bianchi geometries 



Type 


m 


n 2 


n 3 










I 











Type 


m 






II 








+ 


V 











VIo 







+ 


Class B IV 








+ 


VIIo 





+ 


+ 


vi ft 







+ 


VIII 




+ 


+ 


Vllfc 





+ 


+ 


IX 


+ 


+ 


+ 











Table H gives the classification of Bianchi geometries. Note that the invari- 
ance of the Bianchi types under permutations and sign changes of the frame 
elements has been used to simplify the presentation. Here the notation VIo, 
VIIo, Vlfe, VII/j refers to the value of h defined by (|3.2j) . In the list of Bianchi 
types I-IX, the missing type III is the same as VI_i. 

Due to the local homogeneity of the Cauchy surfaces M in a Bianchi 
spacetime, the topologies of the spatially compact Bianchi spacetimes can be 
classified using the classification of compact manifolds admitting Thurston 
geometries. 

The eight Thurston geometries S 3 , E 3 , H 3 , S 2 xR, H 2 x R, Nil, SL(2,R), 
Sol, are the maximal geometric structures on compact 3-manifolds, see |145| 
llfiOj for background. Each compact 3-manifold with a Bianchi (minimal) 
geometry also admits a Thurston (maximal) geometry, and this leads to a 
classification of the topological types of compact 3-manifolds with Bianchi 
geometry, i.e. compact manifolds of the form X/T where X is a complete, 
simply connected 3-manifold with a Bianchi geometry and T is a cocompact 
subgroup of the isometry group of X. It is important to note that V is not 
always a subgroup of the 3-dimensional Bianchi group G. 

The relation between the Bianchi types admitting a compact quotient and 
the Thurston geometries is given by Table EJ For each Bianchi type we give 
only the maximal Thurston geometries corresponding to it, see |74 (I117(I118| 
for further details and references. We make the following remarks 

Remark 3.1. (i) Let (M,g) be a 3-dimensional space form with sec- 
tional curvature k = — 1,0, +1. A spacetime (M,g) with M = 
M x (a, b) and a warped product metric g = —dt 2 + w 2 {t)g, satisfy- 
ing the perfect fluid Einstein equations is called a (local) Friedmann- 
Roberts on- Walker (FRW) spacetime. Specifying the equation of state 
for the matter in the Einstein equations leads to an ODE for w. The 
FRW spacetimes play a central role in the standard model of cosmol- 
ogy. In the vacuum case, only k = —1,0 are possible, and in this 
case, the spatially compact local FRW spacetimes are for k = 0, the 
flat spacetimes covered by T 3 x E 1 , a special case of Bianchi I, and 
for k = —1, the local FRW spacetimes discussed in example \2.,5[ 
which are Bianchi V. 
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Table 2. Bianchi and Thurston geometries 



Bianchi type 


class 


Thurston geometry 


comments 


I 


A 


E a 




II 


A 


Nil 




III = VI_i 


B 


H 2 x E 1 

SL(2,R) 


cf. E.emark IH. IBiiil 


IV 






no compact quotient 


V 


B 




cf. Remark 13.1 Hi vl 


VIo 


A 


Sol 




vift, mo,-i 






no compact quotient 


VIIq 


A 


E 3 




VH h , h + 


B 


H 3 


cf. Remark |3.1||iv| 


VIII 


A 


SL(2,R) 




IX 


A 


S 3 





(ii) The Thurston geometry S 2 x R admits no 3-dimensional group of 
isometries, and hence it does not correspond to a Bianchi geome- 
try, but to a Kantowski-Sachs geometry 105 with symmetry group 
SO (3) x R. The other cases with ^-dimensional symmetry group are 
vacuum LRS Bianchi I, II, III, VIII and IX, see |661 p. 133]. 

(hi) The type of geometry depends on the subgroup T of the isometry 
group used to construct the compactification. The isometry group 
in turn depends on the Bianchi data. Note that 5X(2,R) is both 
maximal and minimal, whereas the isometry group of H 2 x R has 
dimension 4 and is therefore not a minimal (Bianchi) group. 

(iv) The compactifications of Bianchi V and Vllh, h^O are both of the 
type discussed in Examvle \2.h\ thus no anisotropy is allowed in the 
compactification of Bianchi type V and Vllh, h ^ 0. 

(v) The compactification of a Bianchi geometry, introduces new ( moduli ) 
degrees of freedom, in addition to the dynamical degrees of freedom, 
see |117j and references therein for discussion, see also 50 j. The 
resulting picture is complicated and does not appear to have been 
given a definite treatment in the literature. 

In the rest of this section, we will concentrate on class A Bianchi space- 
times. We will also refrain from considering the moduli degrees of freedom 
introduced by the compactification, as it can be argued that these are not 
dynamical. 

Let the expansion tensor 9ij be given by 

@ij — Vj Ui , 

(i.e. 9ij = —kij where kij is the second fundamental form). Decompose 9ij 
as 9ij = o~ij + H8ij where H = 9/3, 9 = g lJ 9ij/3. Then H is the Hubble 
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scalar. We are assuming (J3.1JI and in the vacuum class A case it follows that 
Uij is diagonal |163| p. 41]. 

Under these assumptions, 7^ = 7^,(t) and we may describe the geom- 
etry of (M,g) completely in terms of 7^ or equivalently in terms of the 
3-dimensional commutators 7^ and the expansion tensor given in terms of 

<7jj, 9. 

By a suitable choice of frame, <7jj and riij can be assumed diagonal. Since 
&ij is traceless and diagonal, it can be described in terms of the variables 

and similarly, rtj,- can be represented by 

(ni, 712,713 ) = —/=(7lll, 7122, 7133)- 

The curvature 6^ = Ric(ei,ej) can be written in terms of n^, fry = 2^ nkj— 
n^riij and is therefore diagonal in the chosen frame. Decomposing bij into 
trace-free and trace parts s^,-, fc we find that Sj,- can be represented by two 
variables s+, s_. 

Next we introduce the dimensionless variables (following Hsu and Wain- 
wright |164j . but normalizing with H instead of 9), 

(E + ,E-,N 1 ,N2,N 3 ) = (a + ,a^,n 1 ,n 2 ,n 3 )/H 



Similarly, we set 

(S+,S-) = (s + ,s-)/H 2 

In addition to these choices we also define a new time r by e T = £, t the 
length scale factor, or dt/dr = 1/H. 

Clearly the equations are invariant under permutations (Ej) — > P(Ej), 
(Ni) — > P(Ni). Cyclic permutations of (A^),(Ej) correspond to rotations 
through 2-7r/3 in the E + ,E_ plane. 

We now specialize to the vacuum case. Then (for H 7^ 0) the Einstein 
equations are equivalent to the following system of ODE's for E + , E_, Ni, 
N2, N3 (where the ' denotes derivative w.r.t. the time coordinate r): 



N[ = (q-4£ + )Ni, (3.3a) 

N' 2 = (<? + 2E+ + 2\/3E_)A^ 2 , (3.3b) 

N^ = {q + 2E+ - 2v / 3E_)7V 3) (3.3c) 

E' + = -(2-g)E + -35+, (3.3d) 

E'_ = -(2-g)E_ -35_, (3.3e) 
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where 

q = 2(Z 2 + + Z 2 _), 

S+ = |[(JV 2 - iV 3 ) 2 - Ni(2m -N 2 - N 3 )\, 

S- = -^={N 3 - N 2 )(Nt -N 2 - N 3 ). 

The Hamiltonian constraint (|2.2a|) is in terms of these variables 

£ 2 + S 2 + [Nf + Ni + N$ - 2(N 1 N 2 + N 2 N 3 + N 3 N{)] = 1. (3.4) 

With our conventions, if (r_, t+) is the maximal time interval of existence 
for the solution to (|3.3j) . then r — > r_ corresponds to the direction of a 
singularity. For all non-flat vacuum Bianchi models except Bianchi IX, the 
spacetime undergoes an infinite expansion as r — > r + and is geodesically 
complete in the expanding direction, cf. |134l Theorem 2.1] which covers 
the vacuum case as a special case. For Bianchi IX, on the other hand, 
t — > t _|_ corresponds to H — > 0, as follows from the proof of the closed 
universe recollapse conjecture for Bianchi IX by Lin and Wald |121j . and 
hence to the dimensionless variables becoming ill defined. 

We will now review the basic facts for Bianchi types I, II and IX. 

I: Kasner. The Hamiltonian (|3.4j) constraint reads 

S^ + S 2 =l, (3.5) 

and the induced metric on each time slice is flat. These spacetimes can be 
given metrics of the form 

ds 2 = -dt 2 + t2pt dx i ® dx\ (3.6) 

i 

where ^^Pi = 1, ^2iP 2 = 1 ; the Kasner relations, which correspond to 
the equation (|3.5j) . The H± and pi are related by S + = ^(p 2 + p 3 ) — 1, 

X_ = ^Y-{p 2 —p 3 )- Clearly, the Kasner circle given by (|3.5|) consists of fixed 
points to the system H3.3|) . The points T\,T 2 ,T 3 , with coordinates (—1,0), 
(1/2, ±\/3/2) correspond to flat spacetimes of type I or VIIo (quotients of 
Minkowski space). 

II: By using the permutation symmetry of the equations, we may assume 
N 2 = N 3 = 0. The solution curve which is a subset of the cylinder Yj 2 + +T?_ < 
1 has a past endpoint on the longer arc connecting T 2 and T 3 and future 
endpoint on the shorter arc connecting these points see figure |2j This curve 
realizes the so called Kasner map, cf. 163] §6.4.1]. We define the Bianchi II 
variety as the union of the three spheres in M 5 given by the solutions to the 
constraint (|3.4j) with (N 2 ,N 3 ) = and permutations thereof. The Kasner 
circle is the intersection of the spheres. 
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IX: Mixmaster, characterized by Ni > 0, i = 1, 2, 3. The heuristic picture is 
the following. The projection in the (£+, £_)-plane of a generic orbit in the 
direction r \ r_ moves into the Kasner circle and stays there, undergoing an 
infinite sequence of bounces, which are approximately given by the Kasner 
billiard, cf. figure 0] This picture is supported by numerical studies of the 
full Bianchi IX system, see eg. [30] . 

The Kasner billiard is the dynamical system given by mapping a (non- 
flat) point p on the Kasner circle to the point on the Kasner circle which is 
the end point of the type II orbit starting at p. This map can be described 
as follows. Let B be the nearest corner to p of the triangle shown in figure 
|1J The ray starting at B through p intersects the Kasner circle in a point 
q, which is the image of p under the Kasner map, see also |163l Fig. 6.13]. 
Iterating this construction gives a sequence of points {pi} on the Kasner 
circle, which we may call the Kasner billiard. 

The exceptional orbits which do not exhibit this infinite sequence of 
bounces are the Taub type IX solutions, cf. figure Up to a permuta- 
tion these are given by the conditions 

N 2 = N 3 , £_ = 0. (3.7) 

We call a Bianchi spacetime, satisfying (|3.7|) up to a permutation, a Taub 
spacetime. 

The past limit of the Taub type IX solution is the flat point (—1,0). 
The MVCD of Taub type IX data has a smooth Cauchy horizon, and is ex- 
tendible, the extension being given by the so-called Taub-NUT spacetimes. 
As shown by Chrusciel and Rendall 55, Theorem 1.2], these are the only 
Bianchi IX spacetimes with a smooth Cauchy horizon, which gives a version 
of SCC for this class. See [55., for further details on the status of SCC in the 
locally homogeneous case. Chrusciel and Isenberg proved that the MVCD of 
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Figure 3. A Taub type IX solution, the Kasner circle and 
the triangle for the Kasner billiard in the (£+, S_)-plane are 
shown. 



Taub type IX data has non-isometric maximal vacuum extensions, further 
emphasizing its pathological nature, cf. [53*] . 

From the point of view of the cosmic censorship conjecture, the follow- 
ing theorem appears fundamental. A point which is a past limit point of 
(S + (t), £_ (r), Ni(t), JV 2 (r), N^r)) is called an a limit point. We say that 
the approach to the singularity is oscillatory, if the set of a limit points 
contains at least two points on the Kasner circle, at least one of which is 
disctinct from the special points T\ , T2, T3 , cf . |137| . 

Theorem 3.1 ( |137| I141j ). A vacuum Bianchi spacetime of class A has 
exactly one of the properties 

(1) The Kretschmann scalar n = R a p^sR a ^ & satisfies limsup,-^^ \k\ = 
00 

(2) The MVCD has a smooth Cauchy horizon and the spacetime is a 
Taub spacetime. 

For non-Taub vacuum Bianchi VIII and IX spacetimes, the approach to the 
singularity is oscillatory. 

Remark 3.2. (i) Rendall |137j proved the dichotomy in Theorem I.V. 1\ 
for all Bianchi class A except VIII and IX. These cases as well as the 
oscillatory behavior for type VIII and IX were proved by Ringstrom 
|141j . Curvature blowup for non-Taub Bianchi A models with perfect 
fluid matter, including stiff fluid was proved by by Ringstrom |142j . 
(ii) In class B it is only the exceptional model Bianchi VI-i/g which 
exhibits oscillatory behavior as has been argued by Hewitt et. al |87j . 
(hi) Theorem shows that SCC holds in the class of vacuum Bianchi 
class A spacetimes, also with respect to C 2 extensions. 
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Figure 4. Bianchi IX orbit showing a few bounces. The 
vertical axis is N%. The Kasner circle and the triangle for 
the Kasner billiard in the (£+, E_)-plane are shown. 

(iv) See also Weaver |167j for a related result for Bianchi VIq with a mag- 
netic field. Weaver proved that the singularity in magnetic Bianchi 
VIq is oscillatory and that curvature blows up as one approaches 
the singularity. It should be noted that vacuum Bianchi VIq, on the 
other hand, is non-oscillatory. 

The dynamics of the Bianchi spacetimes has been studied for a long time, 
from the point of view of dynamical systems. In particular, it is believed that 
the Bianchi IX (mixmaster) solution is chaotic in some appropriate sense, 
see for example the paper by Hobill in |163j or the collection jHH] as well as 
the work of Cornish and Levin |57j . for various points of view. However, it 
does not yet appear to be clear which is the appropriate definition of chaos 
to be used, and no rigorous analysis exists for the full Bianchi IX system. 
In the course of the above mentioned work, approximations to the Bianchi 
dynamics have been described and studied, such as the Kasner billiard (cf. 
fig- HI and the discussion above) and the BKL map, cf. |163l §11.2.3]. 

It has long been conjectured that the Bianchi II phase space is the attrac- 
tor for the Bianchi VIII and IX dynamics, see [163 . The phase space for the 
Bianchi II model is given by the Hamiltonian constraint (|3.4j) together with 
the condition N1N2 = N2N3 = N\Ns = 0. This defines a variety consisting 
of the union of three 2-spheres in M 5 . The Kasner circle is the intersection 
of these spheres. The conjecture has been proved in the Bianchi IX case. 

Theorem 3.2 ( 11421 ). The Bianchi II variety is the asymptotic attractor 
for vacuum Bianchi IX. 

This result goes a goes a long way towards proving the chaotic nature of 
the Bianchi IX dynamics. 

The questions of curvature blowup and oscillatory approach to the sin- 
gularity can be studied also in the case of Bianchi class B models. In this 
case we have curvature blowup except for Bianchi III and V. However, in 



THE GLOBAL EXISTENCE PROBLEM IN GENERAL RELATIVITY 



23 




Figure 5. An orbit of the Bianchi VI_i/g system, showing a 
few bounces. Note that in contrast to the Bianchi IX system, 
the dynamics shows a combination of several types of 
bounces. 



class B it is only the so-called exceptional Bianchi VI_i/g which, based on 
numerical work and qualitative analysis by Hewitt et. al [SJj has an oscil- 
latory singularity. The exceptional Bianchi VX_i/g has the same number 
of degrees of freedom as the most general class A models, Bianchi VIII and 
IX. Figure shows an orbit of the Bianchi VX_i/g system. 

4. G 2 

In this section we consider the case when (M,g) is a 3+1-dimensional, 
spatially compact, globally hyperbolic, vacuum spacetime, with a 2-dimensi- 
onal local isometry group G2 with the action of G2 generated by space- 
like Killing fields. By passing to the universal cover, we see that the non- 
degenerate orbits of the G2-action are 2-dimensional homogeneous spaces 
and hence the induced metric on the orbits must have constant curvature. 
The isometry group of the sphere S 2 has no two dimensional subgroup and 
thus the orbits have geometry E 2 or H 2 . 

The special case when the group U(l) x U(l) itself acts on M was con- 
sidered by Gowdy [HUJ|^. In this case, it follows that orbits are compact 
(unless there is an extra Killing field), and M is covered by T 3 , S* 3 or S 1 xS 2 . 
Suppose p is a fix point for the action of U(l) x U(l) on M, then U(l) x U(l) 
acts by isometries on T p M. As the isometry group SO(3) of T p M does not 
have a 2-dimensional subgroup, any degenerate orbit of G must be a closed 
curve. 

Let M be a bundle over S 1 x R with compact 2-dimensional fiber F and 
suppose that the orbits of the G-action on the universal cover of M cover the 
fibers F. If F has geometry E 2 it follows that the Killing fields generating 
the G-action commute, and hence it is natural, following Rendall |136j . to 
use the term local U(l) x U(l) symmetry for this situation. 
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Space-times with local U(l) x U(l) symmetry have also been considered 
by Tanimoto |155j , who discussed the question of in which case a spacetime 
with local U(l) x U(l) symmetry can be considered as a dehomogenization 
of a Bianchi spacetime. It can be seen from the structure of the Bianchi 
groups that the only U(l) x U(l) symmetric spacetimes which have Bianchi 
(or Kantowski-Sachs) limits are covered by T 3 or S 3 , and in the case of S* 3 
(Bianchi IX), it is only the Taub metrics that admit a U(l) x U(l) action 
by isometries. If we consider the case with local U(l) x U(l)-symmetric 
case on the other hand, then all Bianchi models except Bianchi VIII and 
XI, and in case of Bianchi VIII and IX the Taub metrics, can be viewed as 
limits of locally U(l) x U(l)-symmetric models, and therefore these serve 
as dehomogenization of the Bianchi models. 

For simplicity, we concentrate in the rest of this section only the case 
of U(l) x U(l) symmetric spacetimes, and assume that the twist constants 
vanish, i.e. i.e. £1 A £2 A d£i = £1 A £2 A d£2 = 0, where £i,£2 are one 
forms dual to the generators of the U(l) x U(l) action. This is a nontrivial 
restriction only in case M = T 3 , see |81l p. 211]. Such spacetimes are 
known in the literature as Gowdy spacetimes. We further specialize to the 
case with Cauchy surface M = T 3 

Let t,x be coordinates on the 1+1 dimensional Lorentzian orbit space 
(U(l) x U(1))\M and let A(t,8) be the area of the orbit. Gowdy showed 
that in the non-twisted case, with M = T 3 , there are no degenerate orbits, 
i.e. and further, the level sets of A in (U(l) x U(1))\M are space-like. 

We may therefore choose coordinates so that A = 4ir 2 e~ T , and choose the 
metric on the orbit as Ah where h = h(r, x) is a unit determinant metric. 
M is causally incomplete in the direction r f 00, which corresponds to a 
cosmological singularity. 

By construction, the metric on the orbit, h = h(r, x) is a unit determinant 
metric which is constant on each orbit. It therefore represents an element 
of the Teichmiiller space T(T 2 ). The space T{T 2 ) with the Weil-Peterson 
metric is isometric to the hyperbolic plane H 2 . The identification of T(T 2 ) 
with H 2 gives a map u : (U(l) x U(1))\M — * H 2 . This can be realized 
concretely for example by using the model for H 2 with metric 

dP 2 + e 2P dQ 2 (4.1) 

and letting u = (P, Q) with 

P = ln(/m), Q = e ~ P h 12 . 

This is the parametrization that is used in the numerical work of Berger 
and collaborators. Thus the Gowdy (as well as the general U(l) x U(l) 
symmetric) Einstein equations on M = T 3 x R can be viewed as equations 
for the evolution of a loop in H 2 . The velocity of a point u(x,t) in H 2 is 
given by w H 2(x,r) = yj (d T u, d T u). Define the asymptotic velocity t) H 2(x) 
by #jja(a?) = limT-^oo v H 2(x, r), when the limit exists. 
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The spacetime metric may then be written in the form 5 

l^ds 2 = e-V2+r/2 ( _ e -2r dr 2 ^2) 

(4.2) 
+ e~ T [e p dy 2 2 + 2e p Qdy 2 dy 3 + (e p Q 2 + e~ p )dyl] 

where £o is the unit of physical length. Here (r, 9) £ R x S 1 are coordinates 
on the orbit space, y" 4 ,^ = 2,3 are coordinates on the orbit, < y A < 2ir. 
The G invariance implies that we can assume all metric components depend 
on r, 9 only. Let e = dx 2 + (dy 2 ) 2 + (dy 3 ) 2 . The lapse function N satisfies 
N = y / det(g)/ det(e) and hence the time function r is spacetime harmonic, 
cf. subsection 12.31 

Let r\ = —dr 2 + e 2r dx 2 . The Einstein evolution equations take the form 

r] aP {d a pu a + T a bc (u)d a u%u c ) = 0, (4.3) 

where are the Christoffel symbols on H 2 with the metric (|4.1|) . The 
system (|4.3|) is supplemented by a pair of equations for A, which are implied 
by the Einstein constraint equations, and which are used to reconstruct the 
3+1 metric g. 

Equation (|4.3j) is a semilinear hyperbolic system, which resembles the 
wave-map equation, (|2,8|) . Energy estimates or light cone estimates prove 
global existence on (0, oo) x S . The wave operator — d 2 +e~ 2r d 2 degenerates 
as t / oo, which corresponds to a singularity in the 3+1 spacetime (M,g), 
since the area e~ T of the orbit tends to zero. 

The energy E = \ J sl (d T u,d T u) + e~ 2r (d x u, d x u) = Ek + Ey, where 
Efc,Ey are the kinetic and potential energy terms, respectively, satisfies 

d T E = -2E V . 

Therefore E is monotone decreasing, with a rate determined by Ey. This 
shows immediately that there is a sequence of times (r^), lim^oo = 
oo, so that Ey(Tk) — > as k — > oo, which indicates that the scale-free 
variables (e~ T d x P, e~ T d x Q) become insignificant for the dynamics as r — > oo. 
This leads to the idea that the Gowdy system behaves asymptotically as 
a dynamical system in the likewise scale-free variables (d T P,d T Q). This 
heuristic is, as we shall see, supported by numerical work and some rigorous 
results. 

We now briefly discuss some results and open problems for the Gowdy and 
U(l) x U(l) symmetric spacetimes. The following result proves Conjecture 
121 for this class. The first global existence result for Gowdy spacetimes 
with topology T 3 x R was due to Moncrief |125j . who proved that vacuum 
Gowdy spacetimes with the stated topology are globally foliated by level 
sets of the area function, assuming the existence of a compact level set of 
the area function. This assumption was removed by Chrusciel [17] who also 
studied global properties of the area function for Gowdy spacetimes on S s 

5 We have assumed here that the Killing fields are hypersurface orthogonal, see |47| for 
the most general form of the U(l) x U(l) symmetric metric on T 3 x R with vanishing 
twist. 
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and S 2 x S 1 . A class of "nongeneric" metrics still remains to be studied, see 
|47j . The first result concerning global CMC foliations in Gowdy spacetimes 
was proved by Isenberg and Moncrief |99j for the case of vacuum Gowdy 
spacetimes. Recently the following was proved. 

Theorem 4.1 (Andreasson, Rendall, Weaver Nonflat spatially com- 

pact spacetimes with (U(l) xU(l)) symmetry and Vlasov matter, are globally 
foliated by CMC hypersurfaces with mean curvature taking on all values in 
(-oo,0). 

It was conjectured by Belinskii, Khalatnikov and Lifschitz 120 that in 
a generic spacetime with a cosmological singularity, spatial points will de- 
couple as one approaches the singularity, and spatial derivatives become 
insignificant asymptotically. This leads to the idea that asymptotically near 
the singularity, the dynamics of the gravitational field should be explained 
by a family of ODE systems. This is a very rough idea. Yet, as it turns out, 
this principle appears to hold in the cases we have been able to study. The 
idea of BKL has been specialized and reformulated by among others Eard- 
ley, Liang and Sachs jHlj and Isenberg and Moncrief |101j . into the notion 
of asymptotically velocity term dominated (AVTD) singularities. Roughly, 
an AVTD solution approaches asymptotically, at generic spatial points, the 
solution to an ODE, the parameters of which depend on the spatial point. 
In particular, in an AVTD spacetime, locally near a fixed spatial point the 
spacetime approaches a Kasner limit, with parameters depending on the 
spatial point. The asymptotic model for an AVTD spacetime necessarily 
has a non-oscillatory singularity, which excludes generic spacetimes. The 
precise formulation of the BKL conjecture for generic spacetimes is more 
subtle, see [1611 ITS] for discussion of this problem. 

The numerical studies referred to above indicate that for a generic Gowdy 
spacetime, as r /* oo, the velocity u H 2 is eventually forced to satisfy < 
"^H 2 — 1) except at isolated x-values, even if u H 2 > 1 in some subsets of S 1 
initially, and further that u H 2 has a limiting value £>jj 2 (x) for each x as one 
moves toward the singularity. 

The numerical solutions exhibit "spikes" at those x-values, where v h 2 > 1 
asymptotically, cf. figured The very sharp spikes seen (in Q) are coordinate 
effects, corresponding to a part of the solution loop approaching the point 
on the boundary of H 2 sent to infinity by the transformation leading to the 
model with metric Q4.1JI . 

Gowdy spacetimes such that d x Q = are called polarized. Equation 
(|4.3j) then becomes linear. It was proved by Isenberg and Moncrief |101j 
that polarized Gowdy spacetimes are AVTD. Numerical studies by Berger 
et. al., see [2^] and references therein, see also (EHIES], support the idea 
that general Gowdy spacetimes are AVTD. It should be noted that there 
are polarized Gowdy spacetimes with v > 1 up to the singularity, but the 
above mentioned work indicates that this behavior is non-generic. See also 
|48j for work on the case v = 0. 
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P and Q at time T == 40 




1 2 3 4 5 6 



Figure 6. Spikes in P,Q. The very sharp spikes (in Q), 
so-called "false spikes" are coordinate effects. 

The equation (|4.3|) may, essentially under the restriction < v < 1, be 
written as a Fuchsian system 

(tdt + E)U = F{U], 

where t = e~ T , and hence using a singular version of the Cauchy-Kowalew- 
skaya theorem, cf. Kichenassamy and Rendall |1U7| and references therein, 
see also \V3\, AVTD solutions may be constructed given real analytic "data 
on the singularity". The solutions constructed are in terms of the time 
coordinate t, of the form 

P(t, x) = v(x)t + <j)(x) + e~ eT u(x, r) (4.4a) 

Q(r, x) = q(x) + e ~ 2v ^ T [if>(x) + w(t, x)] (4.4b) 

where e > 0, < v < 1, and !«,u^0asr^ do. The Fuchsian method was 
generalized for the Gowdy case by Rendall [I38 a to the C°° case. 

Rendall and Weaver |14()j have constructed families of solutions with 
spikes, starting from a solution with a false spike at xq with k = Vh2(xq) < 1, 
by applying an explicit transformation, a new solution with u H 2(xo) = 1 + k. 
In the new solution £ H 2 has a discontinuity at xo with lim^^^,, z) H 2 (x) = 
1 — k. By iterating this procedure spikes with arbitrarily high velocity can 
be constructed. The group formed by composing the transformations used 
by Rendall and Weaver has been termed the Geroch group [HE] ■ Spikes with 
asymptotic velocity z) H 2 > 2 correspond to higher order zeros of d x Q, and 
are therefore non-generic. Garfinkle and Weaver [JZj have studied the dy- 
namics of generic spikes with initally high velocity. Their analysis, based 
on a combination of heuristic arguments, related to the "method of con- 
sistent potentials", see (2H1 for discussion, and numerical work, shows that 
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the velocity of these spikes is driven down into the interval 1 < v < 2 by a 
sequence of bounces, the qualitative features of which can be explained in 
terms of the relative importance of the terms in the evolution equation. 

Recently conditions on initial data (essentially a small energy condition 
together with restriction on the velocity) have been given by Ringstrom 
|144[ 114=3] , and Chae and Chrusciel [HE] under which the solution is of the 
form (|4.4|) . This work shows that AVTD behavior holds on open and dense 
subsets of S , and makes earlier work jH2j on formal expansions for the 
Gowdy field equations rigorous. The definition of AVTD solution used by 
Chae and Chrusciel [HE1 Eq. (3.8-3.9)] is more general than (|4.4|) used by 
Ringstrom, allowing more general velocity v but with less precise control on 
the lower order terms. 

In view of the above, it is reasonable to make the following conjecture. 

Conjecture 4. Generic vacuum, spatially compact U(l) x U(l) -symmetric 
spacetimes with vanishing twist are AVTD at the singularity, in the comple- 
ment of an at most countable closed subset E of S 1 . The asymptotic velocity 
w H 2(x) exists for all x G S 1 and satisfies < v(x) < 1 for x G S 1 \ E. u H 2 
is continuous on S 1 \E. 

Remark 4.1. (1) This is implicit in Grubisic and Moncrief |82j . 

(2) Chae and Chrusciel [HE] have constructed solutions with an asymp- 
totic velocity which is discontinuous on any closed set F C S , with 
nonempty interior. In particular F may have nonzero measure. If 
the above Conjecture is correct, this behavior is nongeneric. 

The question of cosmic censorship for the Gowdy spacetimes may be stud- 
ied by analyzing the behavior of the Kretschmann scalar k as r — * oo. 
For the class of polarized Gowdy spacetimes, this was done by Chrusciel 
et. al. |S1]. It was proved by Kichenassamy and Rendall jlU7j that for 
generic AVTD spacetimes constructed using the Fuchsian algorithm, the 
Kretschmann scalar k blows up at the singularity and hence generically, 
these spacetimes do not admit extensions, see also the discussion in [52l 
§3-4]. The same behavior was shown to hold also for the class of Gowdy 
spacetimes with asympotic velocity < v < 1 and bounds on the energy 
density, see |143| l3*H] . Further, the Gowdy spikes constructed by Rendall 
and Weaver also have k blowing up, though at a different rate than nearby 
points. Therefore a reasonable approach to the SCC in the class of Gowdy 
spacetimes, is via conjecture H] 

It is relevant to mention here that the AVTD behavior for Gowdy sym- 
metric spacetimes may be broken by the introduction of suitable matter, 
cf. |168j . where numerical evidence for an oscillatory approach to the sin- 
gularity is presented, for a locally U(l) x U(l) symmetric spacetime with 
magnetic field. Similarly, it is expected that scalar field or stiff fluid matter 
changes the oscillatory behavior of general U(l) x U(l) symmetric models 
to AVTD. However, even in the presence of a scalar field, one expects to see 
spikes forming. 
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Orbit of asymptotic Gowdy system 



Orbit of the full Gowdy system 
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Figure 

7. Projection 

in Kasner plane 

of an orbit of the 

asymptotic Gowdy 

system (|4.5H4.6[t . 
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Figure 

8. Projection 

in Kasner plane of 

the orbit of a point 

of the full Gowdy 

system. 



A systematic way of deriving the system of ODE's governing the as- 
ymptotic behavior of the Gowdy spacetimes is to write (|4,3j) in first order 
form, using the scale invariant operators d T ,e~ T d x , and then cancelling the 
terms involving e~ T d x . Following this procedure and working in terms of H- 
normalized orthonormal frame variables adapted to the U(l) x U(l) orbits, 
constructed in a completely analogous way to those discussed in section |3] 
gives, after cancelling the terms corresponding to e~ T d x , the system 



p s+ 


= (g-2)(l + S + ) 


(4.5a) 




= (q - 2) S_ - 2V3N 2 + 2v / 3£ 2 < 


(4.5b) 


V N X 


= (g + 2£ + )A x 


(4.5c) 




= (q - 2 - Ex - 2V3N X iV_ 


(4.5d) 


V N_ 


= (q+ 2S+ + 2v / 3S_) AL + 2V3 S x N x . 


(4.5e) 



where T>o 
equation 



(1 + X + )<9 r , q = 2(Ei + T?_ + E^.), subject to the constraint 



1 = Si + S 2 + £^ + Nl + N 2 



(4.6) 



See |13j for discussion of the system ()4.5H4.6|) . This system contains the 
spatially homogenous Bianchi class A models I-VIIo, but represented in 
terms of a non-Fermi-propagated frame. 

The asymptotic behavior of the system (|4.5H4.6j) in the direction towards 
the singularity is easily analyzed, and one shows that all solutions of the 
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system have limit points on the Kasner half-circle + X 2 . = 1. Similar to 
the standard Bianchi case, the Kasner circle consists of fixed points. For the 
Bianchi types II-VIIo, only part of the Kasner circle consists of stable fixed 
points. In the present case, only the arc with £_ < 0, S + > \ is stable, and 
a generic orbit of (|4.5|) ends on the stable arc. The stable arc corresponds 
to orbits with asymptotic velocity < -0 < 1. See figure for an example of 
an orbit of the system (|4.5H4.6j) . 

According to the BKL proposal, a generic solution u(t, x) of the full 
Gowdy system will have the property that it will asymptotically approach a 
solution of the Gowdy AVTD system (j4.5H4.6j) . That this is indeed the case 
is indicated by figure |SJ which shows the projection in the Kasner plane of 
the orbit of one point for the full Gowdy system. We see in the figure the 
orbit of a point near a spike point. The points marked on the Kasner circle 
are (from the left) the points with velocity v = do, 2, 1, 0. The stable arc on 
the Kasner circle is plotted with a heavy line. It should be noted that the 
last part of the orbit which ends on the stable Kasner arc is approximately 
a Bianchi II orbit, with initial velocity in the interval (1,2). This is the 
generic behavior, also close to spike points, see |13| for more details. 

A class of AVTD polarized Gowdy spacetimes with non-vanishing twist 
has been constructed by Isenberg and Kichenassamy 96 using the Fuchsian 
algorithm. On the other hand for general non-polarised twisted U(l) x U(l) 
symmetric spacetimes, work of Berger et al. |31j supports the conjecture that 
these spacetimes show oscillating behavior as one approaches the singularity. 
This idea is also supported by the fact that the only Bianchi model with 
U(l) X U(l) symmetry and non-vanishing twist constants is the exceptional 
Bianchi VI_i/g which has been shown heuristically by Hewitt et al. [HI] to 
have oscillatory behavior at the singularity. 

5. U(l) 

The U(l) symmetric vacuum 3+1 Einstein equations is an important case 
which is of intermediate difficulty between the full 3+1 Einstein equations 
and the highly symmetric Gowdy equations. In the presence of a hyper- 
surface orthogonal space-like Killing field, the Einstein equations reduce to 
2 + 1 gravity coupled to wave map matter, the field equations and their re- 
ductions have been derived in |127|I126|I128| . see also [31] ■ Choquet-Bruhat 
and Moncrief have proved that the Cauchy problem for the U(l) problem is 
well-posed in H 2 , see @2j Ell- 
in the above quoted papers, the spacetime is assumed to be a U(l) bundle 
over a spatially compact 2 + 1 spacetime. The case of local U(l) symmetry 
does not appear to have been considered in connection with the Einstein 
equations, see however |131| for information about 3-manifolds with local 
U(l) action. 

It is also possible to study the case when the reduced space is asymptot- 
ically flat. This case has been considered in work by Ashtekar and others, 
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|17| 115] , see also |16] . In these papers an analogue of the ADM mass at spa- 
tial infinity is introduced. It is proved that it is nonnegative and bounded 
from above. It is interesting to study the consequences of the presence of this 
conserved quantity for the 2+1 dimensional Einstein-matter system given 
by the U(l) problem, as one expects that it gives a stronger bound on the 
fields than in the 3+1 case. This appears to be a natural setting for a small 
data version of the U(l) problem. 

In the following, we will consider the spatially compact case. Let (M, g) 
be a 3+1 dimensional spacetime, assume M = B x R, with ir : B — > S a 
principal U(l) bundle, £ a compact surface. Further assume the group U(l) 
acts by isometries on (M,g) with the action generated by the Killing field 
J, which we assume to be space-like, (J, J) > 0. 

Let the function A on £ x R be defined by ir*\ = | log ((J, J)), and let 
9 = e~ 2n x J. Then we can write the spacetime metric g in the form 

2 = 7r *( e - 2A 5 ) + e 2 -* A 6>«>6>, 

where g is a Lorentzian metric on R x S. 

Introduce a frame 6 e a , a = 0, 1, 2, 3 with = e~ x J, and let the indices 
a, b, c, • • • = 0, 1, 2. We may without loss of generality assume that [J, e a ] = 0. 
We have 

dJ a p = O a/3 + 2(e a (\)Jp - e /3 (A) J Q ), 

where ® a pJ^ = 0. It follows that @ a p = 7r*(e 2X F a p) where F a g is a 2-form 
on R x S. If Fq,^ = 0, then J is hypersurface orthogonal. 

To avoid cluttering up the notation we will in the following make no 
distinction between fields on the orbit space Mx5 and their pullbacks by 

7T. 

The components of the 3+1 Ricci tensor are 

Rab = Rab - 2V a AV 6 A + V c V c Xg ab - ^e 4X F ac F a c , 
R a3 = ^e- A V c (e 4A F a c ), 

R33 = e 2X [-g ab V a V b X + -^ x F ab F ab ] (5.1) 

Let the one— form E on R x S be given by E = — * g (e 4A i ? ). One of the 
Einstein equations (|5.1]> implies that dE = 0. 

Now the Einstein vacuum equations R a p = imply the system 

Rab = ^(4V a AV fe A + e~ AX E a E b ), (5.2a) 

V a V a A + X -e~ AX E a E h g ah = 0, (5.2b) 
V a V a a; - 4V a XE a = 0. (5.2c) 

6 In |41) . a frame dx a ,0 3 is used, where 8% — e~ 2X J a = e _A e^ . Furthermore, their 

F a f3 = dOa/3 — e~ 2X Q a/ 3. 
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When E a = V a w, we recognise (|5.2jl as the 2+1 dimensional Einstein 
equations with wave map matter, for the wave map with components (A,w) 
with target hyperbolic space H 2 with the constant curvature metric 



See subsection 12.41 for some discussion of wave map equations. 

We now specialise even further to the polarized case E a = 0. This corre- 
sponds to assuming that the bundle it : B —* £ is trivial, and that the vector 
field J is hypersurface orthogonal. Then to is constant and the equations 
(|5.2j) become 



which is precisely the 2+1 Einstein equations coupled to a massless scalar 
field. 

In the 2-dimensional case, the operator k — > divfc is elliptic on symmetric 
2-tensors with vanishing trace. Further, by the uniformization theorem, a 
compact 2-dimensional Riemannian manifold (£, h) is conformal to (£, [h]) 
where [h] is a representative of the conformal class of g, i.e. a constant cur- 
vature metric. Therefore working in spatial harmonic gauge with respect to 
(£, [h]) (i.e. conformal spatial harmonic gauge) and CMC time gauge, the 
constraint equations form an elliptic system for (hij, kij). We get a represen- 
tation of (hij, kij) in terms of ([h], 4>, fc TT , Y) where [h] is the conformal class 
of h, corresponding to a point in Teichmiiller space, <j> a conformal factor de- 
termined by the Hamiltonian constraint equation (a nonlinear elliptic system 
for (/>), fc TT a trace-free, divergence free 2-tensor on (£, [h]), corresponding 
to a quadratic differential, and finally, Y is a vector field determined from 
the momentum constraint equation. Note that ([/i],fc TT ) represents a point 
in T*T(S), the cotangent bundle of the Teichmiiller space of E. 

Due to the ellipticity of the constraint equations in the 2+1 dimensi- 
onal case, in the gauge as described above, it is possible to eliminate the 
Einstein equation (|5.2a|) from the system Q5.2JI and instead solve the elliptic- 
hyperbolic system consisting of the hyperbolic system (|5,2bl5.2cj) . coupled 
to a nonlinear system which determines components of the spacetime metric 
g in terms of the data via the constraint and gauge fixing equations, as well 
as an ODE system determining the evolution of the Teichmuller degrees of 
freedom of the metric on E. In the special case £ = S 2 , Teichmuller space is 
a point, and further H 1 (S 2 ) = 0, which means that E = duj. Therefore ig- 
noring problems with gauge fixing, the reduced U(l) system would consist of 
the wave equation coupled to the system which determines the components 
of g. 

A special case of the polarized U(l) equations is given by setting A = 
constant. Then the field equations (|5.^j) are just 




Rab = V a AV 6 A 
V a V a A = 



(5.3a) 
(5.3b) 



Rab = 0, 
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the 2+1 dimensional vacuum equations. In this case, the spacetime is a 
3-dimensional Lorentzian space-form. The dynamics of 2+1 dimensional 
vacuum gravity has been studied by Andersson, Moncrief and Tromba 
who proved global existence in CMC time for 2 + 1 dimensional vacuum 
spacetimes, with cosmological constant, containing at least one CMC hyper- 
surface. This proves Conjecture El for the class of 2+1 dimensional vacuum 
spacetimes. 

Andersson and Rendall ^2j showed that the 3+1 dimensional Einstein 
scalar field system can be formulated as a Fuchsian system. Using similar 
techniques, Damour et al. |611 §4] have shown that the polarized U(l) equa- 
tions may be formulated as a Fuchsian system, and therefore AVTD solutions 
may be constructed using a singular version of the Cauchy-Kowalewskaya 
theorem, as was done by Kichenassamy and Rendall jlU7j for the Gowdy 
case. See also Isenberg and Moncrief |1(J2| which deals with the "half- 
polarized" U(l) system in addition to the polarized case. This supports the 
following conjecture. 

Conjecture 5. Generic polarized U(l) spacetimes are AVTD at the singu- 
larity. 

Remark 5.1. Conjecture\^was essentially stated by Grubisic and Moncrief 
|83j . This is supported by numerical work of Berger and Moncrief |33j . 
Similarly to the case of Gowdy spacetimes, see Conjecture^ one expects that 
" 'spiky features" will form at the singularity, and hence the AVTD property 
of the U(l) symmetric spacetimes should be understood to be generic also in 
the spatial sense. 

It seems reasonable to expect that polarized U(l) spacetimes which are 
AVTD at the singularity have a strong curvature singularity generically, and 
therefore that proving the AVTD conjecture for polarized U(l) would be a 
big step towards proving SCC for this class. 

In contrast to polarized U(l), the generic U(l) spacetimes have sufficiently 
many degrees of freedom that one expects them to satisfy the BKL picture 
of an oscillatory approach to the singularity, as is also expected in the fully 
3+1 dimensional case, cf. the remarks in section [7| This is supported by 
the numerical evidence so far, see |32| . 

In the expanding direction, on the other hand, a small data semi-global 
existence result holds [421 139j . similar to the one for the full 3+1 dimensional 
case discussed below in section |HJ In this case, the notion of "small data" is 
taken to be data close to data for a background spacetime with space-like 
slice M = X x S 1 where £ is a Riemann surface of genus > 2 with metric 
a of scalar curvature —1, and M has the product metric. In this case, the 
background spacetimes are products of flat spacetimes as in Example 12.51 
with n = 2, with the circle (such spacetimes are of type Bianchi III). The 
proof uses energy estimates for a second order energy, which controls the 
H 2 x H 1 norm of the data, for small data. The energy expression used is a 

2 

combination of several terms. The first order part is ^Vol 5 (S) + 27rx(S), 
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which via the constraint equations may be related to the wave map energy. 
The second order term is defined in terms of the L 2 norm of the Laplacian of 
the wave map field, and the gradient of its time derivative. In order to use 
this method, one of the technical problems one has to deal with is controlling 
the spectral gap of the Laplacian. This is done by keeping control of the 
Teichmuller parameters during the evolution. For technical reasons, the 
assumption is made that the spectral gap Ao of the Laplacian for the inital 
metric satisfies Ao > 1/8. The proof shows that the Teichmuller parameter 
a(t) converges to a point in the interior of Teichmuller space and that the 
spacetime geometry converges in the expanding direction to one of the model 
spacetimes described above. 

6. 3+1 

In the 3+1 dimensional case with no symmetries, the only known facts on 
the global properties of spacetimes are Lorentzian geometry results such as 
the Hawking-Penrose singularity theorems |132j and the Lorentzian splitting 
theorem of Galloway |75j , see |24| for a survey of Lorentzian geometry. 

Here we are interested in results relevant to the SCC, Conjecture ^ and 
the CMC conjecture, Conjecture i.e. results about the global behavior 
of solutions to the Cauchy problem for the evolution Einstein equations, 
in some suitable gauge. The Cauchy problem for the Einstein evolution 
equations has been discussed in section 12.21 

With this limitation there are essentially 3 types of results known and 
all of these are small data results. The results are those of Friedrich on 
the "hyperboloidal Cauchy problem" , of Christodoulou and Klainerman on 
the nonlinear stability of Minkowski space (generalized by Klainerman and 
Nicolo jlllj to exterior domains) and the work of Andersson and Moncrief 
[Sj) on global existence to the future for data close to the data for certain 
spatially compact flat k = — 1 (local) FRW spacetimes, again a nonlinear 
stability result. We will briefly discuss the main features of these results. 

The causal structure of a Lorentz space is a conformal invariant. This 
leads to the notion that the asymptotic behavior of spacetimes can be studied 
using conformal compactifications or blowup. The notion of isolated system 
in general relativity has been formalized by Penrose in terms of regularity 
properties of the boundary of a conformally related spacetime (M,g), with 
null boundary J, such that M is a completion of M, and g = <$>*(Q~ 2 g), 
where $7 6 C°°(M) is a conformal factor, $ : M — ► M is a diffeomorphism 
of M to the interior of M. Given assumptions on the geometry of (M, g) at 
3, Penrose proved using the Bianchi identities that the components of the 
Weyl tensor of (M,g) decay at physically reasonable rates. 

Friedrich derived a first order symmetric hyperbolic system from the Ein- 
stein equations, the "regular conformal field equations". This system in- 
cludes among its unknowns, components of the Weyl tensor, the conformal 
factor Q and quantities derived from the conformally rescaled metric g. This 
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system has the property that under the Penrose regularity conditions at J, 
the solution can be extended across J. The fact that the regular confor- 
mal field equations gives a well posed evolution equation in the conformally 
compactified picture enabled Friedrich to prove small data global existence 
results by using the stability theorem for quasi-linear hyperbolic equations. 

In |69[ Theorem 3.5] Friedrich proved global existence to the future for 
data (M,g,k) close to the standard data on a hyperboloid in Minkowski 
space, satisfying asymptotic regularity conditions compatible with a Penrose 
type compactification (the hyperboloidal initial value problem). This was 
later generalized to Maxwell and Yang-Mills matter in [70]. Initial data for 
the hyperboloidal initial value problem were first constructed by Andersson, 
Chrusciel and Friedrich jBj, see also H], E3- 

The result of Friedrich is a semi-global existence result, in the sense that 
the maximal vacuum Cauchy development D(M) of the data (M,g,k) is 
proved to be geodesically complete and therefore inextendible to the future, 
but not to the past. In fact typically D{M) will be extendible to the past 
and (M, g, k) may be thought of as a partial Cauchy surface in a larger 
maximal globally hyperbolic spacetime (M,g), cf. figure E3 In view of the 
inextendibility to the future of D(M), the result of Friedrich may be viewed 
as supporting the cosmic censorship conjecture. In the case of the Einstein 
equations with positive cosmological constant, the method of Friedrich yields 
a global existence result for data close to the standard data on M = S 3 in 
deSitter space, cf. jSH Theorem 3.3]. 




Figure 9. The semi-global existence theorem of H. Friedrich 

The first true small data global existence result for the vacuum Einstein 
equations, was proved by D. Christodoulou and S. Klainerman They 
proved that for data (M, g, k) sufficiently close to standard data on a hy- 
perplane in Minkowski space, with appropriate decay at spatial infinity, the 
MVCD (M, g) is geodesically complete and therefore inextendible, cf. figure 
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1101 The Christodoulou-Klainerman global existence theorem therefore sup- 
ports the cosmic censorship conjecture. Klainerman and Nicolo |110| llllj 




Figure 10. The global existence theorem of D. Christodou- 
lou and S. Klainerman 

have proved a global existence result for data (M, g, k) which are close to 
standard flat Minkowski space data on an exterior region M \ K, where 
K C M is compact so that M \ K = M 3 \ Ball, and with weaker asymptotic 
conditions compared to the Christodoulou-Klainerman theorem. The result 
of Klainerman and Nicolo states that for a vacuum data set (M,g, k) which 
is sufficiently close to the standard flat Minkowski data on M\K, the outgo- 
ing null geodesies in the causal exterior region D{M \ K) are complete and 
D{M \ K) is covered by a double null foliation, with precise control over 
the asymptotics. It should be noted that the KN theorem therefore cov- 
ers a more general class of spacetimes than the Christodoulou-Klainerman 
theorem and is not strictly a small data result. In particular, the maximal 
Cauchy vacuum development of (M , g, k) may be singular for data satisfying 
the assumptions of the KN theorem, cf. figure 1111 If the smallness assump- 
tion is extended also to the interior region, the Christodoulou-Klainerman 
theorem is recovered. The Einstein equations are quadratic in first deriva- 
tives, and therefore in 3+1 dimensions, one needs something like a null 
condition in order to get sufficient decay for a global existence argument. 
The Einstein equations are not known to satisfy a null condition. However, 
by a detailed construction of approximate Killing fields and approximate 
conformal Killing fields, Christodoulou and Klainerman are able to control 
the behavior of components of Bel-Robinson tensors constructed from the 
Weyl tensor and its derivatives, and close a bootstrap argument which gives 
global existence for sufficiently small data. As part of this argument, it is 
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Figure 11. The exterior global existence theorem of S. 
Klainerman and F. Nicolo. The figure illustrates a situation 
which is covered by this theorem, with a singularity forming 
due to a strong gravitational field in the interior region I(K) 
while the exterior region D(M \ K) has complete outgoing 
null geodesies reaching 3. 

necessary to get detailed control over the asymptotic behavior of light cones. 
Christodoulou and Klainerman also study the asymptotic behavior of com- 
ponents of the Weyl tensor and are able to prove that some, but not all of 
these have the decay implied by the Penrose conditions on J. 

The question about existence of vacuum asymptotically Minkowskian 
spacetimes with regular conformal completion, in the sense of Penrose, has 
been open for a long time. Recently, Delay and Chrusciel were able to 
use a gluing construction of Corvino and Schoen [25 to construct vacuum 
spacetimes which are exactly Schwarzschild near spatial infinity, and which 
are close enough to Minkowski space, so that the global existence result of 
Friedrich applies to show existence of a global J. Further, Klainerman and 
Nicolo |112j have been able to show that a class of asympotically Euclidean 
initial data with suitable decay at infinity has a Cauchy development where 
the null components of the Weyl tensor have decay corresponding to the 
peeling conditions of Penrose. 

The work in [Hj and [3] shows that generic hyperboloidal data do not 
satisfy the required regularity at the conformal boundary, which therefore 
may be viewed as an indication that Penrose regularity at J is non-generic. 
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It is still an open question what general conditions on initial data on an 
asymptotically flat Cauchy surface give a Cauchy development with regular 
conformal completion. H. Friedrich has been developing a programme which 
approaches this problem by analyzing the conformal structure at spatial 
infinite in detail, see [201 for references on this, and see also the paper |162j 
by Valiente-Kroon, which points out some new obstructions to regularity. 

In contrast to the results of Friedrich and Christodoulou-Klainerman- 
Nicolo, the semi-global existence result of Andersson and Moncrief [2] deals 
with spatially compact vacuum spacetimes. 

Let (M, 7) be a compact hyperbolic 3-manifold with metric 7 of sectional 
curvature —1. Then the spacetime M = M x M with metric 7 = — dt 2 + i 2 7 
is a k = — 1 (local) FRW spacetime. A flat spacetime metric on M defines 
a geometric structure with group SO(3, 1) x M 4 . The moduli space of flat 
spacetimes with topology M x R has the same dimension as the moduli 
space of flat conformal structures on M, see 0. We say that M is rigid 
if the dimension of the moduli space of flat conformal structures on M is 
zero. This is equivalent to the condition that there are no nonzero traceless 
Codazzi tensors on (M, 7). Kapovich |106J has constructed examples of rigid 
compact hyperbolic 3-manifolds. 

The global existence theorem proved in [5] states that if M is a rigid ma- 
nifold of hyperbolic type, then for a vacuum data set (M, g, k), sufficiently 
close to the standard data in a spatially compact n = — 1 (local) FRW space- 
time, the MVCD (M, g) is causally geodesically complete in the expanding 
direction. It is a consequence of the singularity theorem of Hawking and 
Penrose that (M, g) is singular, i.e. geodesies in the collapsing direction are 
incomplete, cf. figure WH 




Figure 12. The semi-global existence theorem of L. Ander- 
sson and V. Moncrief. 
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The proof uses local well posedness and a continuation principle for the 
Einstein equations in constant mean curvature and spatial harmonic coor- 
dinates (CMCSH) gauge, cf. subsection 12.31 together with energy estimates 
for a first order Bel-Robinson type energy, defined in terms of the Weyl 
tensor and its covariant derivative w.r.t. T. Let W a b c d be the Weyl tensor 
and let *W a b c d be its dual. In a vacuum spacetime W a b c d = Rabcd- Define 
Eij = WiTjT and B^ = *WiTjT- Then E{j and Bij are the electric and 
magnetic parts of the Weyl tensor. These fields are traceless and symmetric 
and satisfy a Maxwell-like system. A computation shows that the zeroth 
order Bel-Robinson energy 

Q = J\E\ 2 + \B\ 2 
(which measures second derivatives of g in L 2 ), satisfies the equation 
dt = -3 / N[(E x E) ■ k + (B x B) ■ k - -(\E\ 2 + \B\ 2 )trk 

JMt 3 

- 2N- 1 V i N{E A B)i}dnM t 
where for symmetric traceless two-tensors A,B, 

2 

(A x B\j = B m j + Aj B m i — —A B mn gij 

and 

{A/\B) i = e i mn A*B kn , 

Combining Q with the a term defined analogously in terms of V^W^cd and 
extracting a scale invariant quantity, leads to an energy £, which bounds 
the Cauchy data in x H 2 for small data. This energy is shown to satisfy 
an energy inequality of the form 

d t t < -28, - C£ 3/2 

Using this one may show that if £ is small initially, then £ will decrease, 
as long as local existence holds. A continuation principle for the Einstein 
equations in H s x H 2 jlUj completes the proof of global existence for small 
data. 

7. Concluding remarks 

The results and numerical studies so far can be argued to fit with the 
so called BKL picture of cosmological singularities, which states that the 
generic singularity should be space-like, local and oscillatory, see j2H] for 
a review. Roughly, one expects that the "points on the singularity" are 
causally separated and that locally at the singularity, the dynamics under- 
goes a chaotic sequence of curvature driven "bounces" interspersed with 
relatively uneventful "coasting" epochs, with Kasner like dynamics (i.e. the 
metric is, locally in space, approximately of the form given in (|3.6|1 ). The 
locality at the singularity, is easily checked for the nonvacuum Friedman- 
Robertson- Walker models of the standard model of cosmology, and is the 
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cause of the so-called "horizon problem" in cosmology. Locality at the sin- 
gularity can be proved for Gowdy and polarized U(l). 

The picture with bounces and coasting is to a large extent inspired by the 
Mixmaster model (Bianchi IX, cf. section |2J) which is known to have an os- 
cillatory singularity, cf. Theorem l3.il and the above scenario may therefore 
be described by the slogan "generic singularities are spacelike and oscilla- 
tory" . By an abuse of language we may say that the oscillatory behavior is 
of "mixmaster type". 

An approach to understanding the asymptotic dynamics near singular- 
ities is to attempt to extract, in a systematic way, a dynamical system 
governing the asymptotic behavior of the Einstein equations. A promis- 
ing approach is based on an analysis of the evolution equations written in 
Hubble-normalized scale-free orthonormal frame variables along the lines 
discussed for the Gowdy system in section 0J The general situation has 
been discussed in |161j . In this setting the Einstein evolution equations ap- 
pear in first order form, and the natural version of the BKL conjecture is 
that d±U —* in the direction of the singularity, where d\ is a scalefree first 
order operator analogous to the operator e~ T d x for the Gowdy system, and 
U denotes the scale free frame variables. 

If this picture is correct, the asymptotic behavior near the singularity can 
be understood by analyzing the asymptotic behavior of orbits of the system 
of ODE's that one gets by cancelling the terms containing d\ in the first 
order form of the Einstein evolution equations. The possibilities that arise 
have been analyzed in It is now important to understand whether 

or not this version of the BKL conjecture is correct by doing numerical 
experiments. No complete results in this direction have been obtained so 
far. 

The mixmaster type behavior may be prevented by the presence of certain 
types of matter such as a scalar field (stiff fluid). In fact work by Anders- 
son and Rendall [12|, shows that the 3+1 dimensional Einstein-scalar field 
equations may be formulated as a Fuchsian system near the singularity, and 
AVTD solutions may be constructed using a singular version of the Cauchy- 
Kowalevskaya theorem. For the polarized U(l) case, the reduced system is 
2+1 dimensional Einstein-scalar field equations. The Fuchsian method has 
been applied to this special case by Isenberg and Moncrief |102| . Damour et 
al. [HJ have given a rather complete study of the application of the Fuchsian 
method to spacetimes of general dimension with matter. 

It appears likely that some aspect of the picture sketched above will be 
relevant for the final analysis of the large data, global behavior of vacuum 
spacetimes. 

Appendix A. Basic Causality concepts 

Here we introduce the basic causality concepts, see |165[ Chapter 8] . |24| 
lEUESj for details. 
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A compact manifold M admits a Lorentz metric if and only if its Euler 
characteristic x(M) vanishes, while all noncompact manifolds admit Lorentz 
metrics. A vector V G TM is called space-like, null or time-like if g(V, V) > 
0, = or < respectively. V is called causal if it is null or time-like. A 
C 1 curve cCMis called time-like (causal) if c is time-like (causal). This 
extends naturally to continuous curves. A hypersurface Mc Mis called 
acausal if no causal curve meets M more than once, and is called space-like 
if the normal of M is time-like. 

Given an acausal closed subset S C M, the domain of dependence 
D(S) is the set of points p £ M such that any inextendible causal curve 
containing p must intersect S. If S C M is a space-like hypersurface and 
D(S) = M, S is called a Cauchy surface in M. If M has a Cauchy surface 
M is called globally hyperbolic. A globally hyperbolic spacetime has a 
global time function, i.e. a function t on M so that Vt is time-like, with 
the level sets of t being Cauchy surfaces. In particular, M = M x R where 
M is any Cauchy surface. 

M is called time oriented if there is a global time-like vector field on M. 
A globally hyperbolic spacetime is time oriented by the above, and hence 
it makes sense to talk about the future and past domains of dependence 
D + (M) and D~(M) of a Cauchy surface M. In a time oriented space time, 
the chronological future I + (S) of S C M is the set of all points reached by 
future directed time-like curves, of nonzero length, starting on S. The time- 
like past I~(S) is defined analogously. The causal future and past J + (S) and 
J~{S) is defined analogously to 1 (S) with causal curve replacing time-like 
curve (the causal curve is allowed to be trivial). 

The variational problem for geodesies (w.r.t. Lorentzian length) is well 
behaved precisely when (M,g) is globally hyperbolic. Let C(p, q) be the set 
of continuous causal curves between p,q £ M. Then C(p, q) is compact w.r.t. 
uniform convergence for all p,q £ M, if (M,g) is globally hyperbolic. This 
makes global hyperbolicity a natural assumption in Lorentzian geometry, 
just as completeness is a natural assumption in Riemannian geometry. 

Global hyperbolicity may also be characterized by strong causality and 
compactness of C(p, q) for all p,q £ M, or compactness of J + (p) n J~(q) for 
all p, q £ M, see |165| Chapter 8] for details. 

Given a spacetime (M, g) and a space-like hypersurface M C M, the fu- 
ture Cauchy horizon H + (M) is given by H + (M) = D+(M) \ I~[D+(M)], 
with the past Cauchy horizon H~ (M) defined analogously. The set H{M) = 
H + (M) U H~(M) is called the Cauchy horizon of M. It can be proved that 
H{M) = dD(M), the boundary of D(M). 
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